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NONSTANDARD  I)  I  STK  1 IH’T  1 O'.'S  IN 
TWl  -STACK  !.KaS  !'  SOCARKS 

Mark  Mareucci  and  l).K.  lenson 
Virginia  I’o  1  y  tecl.n  ic  Institute 
and  Stat  e  I'n  i  vers  i  t  v 

0.  ABSTRACT 

In  estimating  the  coefficient  of  an  endogenous  variable  in  a  single 
equation  of  a  system  of  Linear  equations,  Anderson  and  Sawn  (1973)  ex¬ 
pressed  the  distribution  of  the  two-stage  least -squares  (TSI.S)  estimator 
ns  o  doubly  nonconi  ra  I  K  distribution.  l\'e  relax  their  assumption  of  in¬ 
dependent  Gaussian  errors,  taking  instead  a  scale  mixture  of  spherical 
Caussian  laws  in  a  class  containing  the  spherical  stable  distributions. 
The  resulting  distribution  of  the  TSI.S  estimator  is  a  mixture  of  double 
noncentral  F  distributions  nixed  over  the  noncentralitv  parameters,  and 
for  suitable  mixtures  the  normal-theory  distribution  is  robust.  Computa¬ 
tions  reported  for  contaminated  spherical  Caussian  and  spherical  Cauchv 
errors  are  compared  with  tilt  standard  case. 

1.  INTRODUCTION 

Various  procedures  have  been  advocated  for  estimating  the  coelficent 
in  simultaneous  systems  of  Linear  equations.  Arguments  supporting  these 
procedures  often  appeal  to  asymptotic  properties  of  the  estimators;  i  •: 
practice,  however,  their  small-sample  proper  ies  arc  o if  interest  as  wo  1  1  . 
A  recent  survey  of  applicable  small-sample  distribution  theory  was  under¬ 
taken  by  Mariano  (1980).  One  technique  in  wide  usage  is  two-stage 
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least-squares  (TS1.S)  estimation  under  t  lie  s  l  and.i  id  .K:uim;'l  inn  that  errors 
of  the  model  are  independent,  ident  i  cn  I  I  v  d  i  s  t  r  ilxit  ed  (iid)  Caussian 
variables.  Under  these  assumptions  Andersen  and  Sawa  ( 1 9 7  i ,  1  '< 79 )  derived 
and  later  tabulated  the  exact  distribution  el  the  ISI.S  estimator  for  the 
coefficient  of  an  endogenous  variable  in  a  single  equation  ■-!  a  linear 
system.  Raj  (1980)  studied  effects  of  mm-faussiaii  errors  on  various 
estimators,  including  the  TSI.S  estimator,  using  Monte  Carlo  methods.  fo 
date,  however,  analytical  investigations  of  the  distributions  >f  ISIS 
estimators  appear  to  have  been  confined  to  the  standard  case  ot  i id 
Caussian  errors. 

In  this  paper  we  study  the  exact  distribution  of  the  TSI.S  estimator 
lor  the  coefficient  of  an  endogenous  variable  in  a  single  equation  under 
nonstandard  assumptions,  replacing  normality  with  an  error  distribution 
in  the  class  of  scale  mixtures  of  spherical  Caussian  laws  (cf.  Kelker 
(.1970),  lor  example).  Our  reasons  for  considering  this  class  follow. 

(i)  beginning  with  Mandelbrot  (I960, |9hj, 1969) ,  infinite-variance 
distributions  including  the  stable  laws  have  been  used  to  model  economic 
variables  such  as  stock  prices;  cf.  also  Kama  (  |9(>r>)  ,  McCullough  (  1 9  7  S )  , 
Kress  (1975),  and  Samuel. son  (1976).  The  class  considered  here  contains 
the  spherical  stable  laws  and  other  heavy-tailed  distributions  whose 
marginals  figure  prominently  in  the  work  cited. 

(ii)  I’lie  d  is  t  r  i  hut  ions  of  Anderson  and  Sawa  (1973,1979)  give  approxi¬ 
mations  in  the  case  of  non-Caussian  errors  attracted  to  the  Caussian  law 
under  conditions  for  central  limit  theory.  I  he  present  sludy  yields 
approximations  in  the  case  of  errors  not  having  second  moments  whose  dis¬ 
tributions  are  attracted  to  a  spherical  stable  law. 

(iii)  Our  approach  admits  a  Bayesian  view  in  which  the  distributions 
studied  are  posterior  distributions  of  the  TSI.S  estimators,  the  errors 


being  cond  i  l  iona  1 1  y  iid  (Iaussian  wif.li  .1  scale  pa  rami'  1  c  r  i!c  1 c  rri  i  ncd  bv  the 
i\m  1  i  /at  i on  of  a  random .  envi  ronmcn t  . 

(  Lv)  The  ordinary  1 easl-s<|ua res  oslimators  in  a  general  I  inoar  modal 
arc  known  to  retain  certain  optimal  propert  i  es  under  spherical  errors  with 
and  without  moments  (rf.  Jensen  (1979)). 

(v)  For  critical  linear  systems  statistical  control  charts  ni.iv  he 
devised  for  monitoring  the  stationaritv  of  a  st  ructurai  relation  based  on 
(iaussian  error  processes  anil  the  work  of  Anderson  and  Sawa  (  1 11 7  3 )  .  The 
present  study  lays  the  foundations  for  developing  such  procedures  in  the 
cast.'  of  spherically  invariant  processes,  as  such  processes  mav  he  representt 
as  stale  mixtures  of  spherical  (Iaussian  processes  (cf.  Hartman  and  Wintner 
(  I  9  ■'»()) )  . 

Our  program  of  study  follows.  In  Section  we  demonstrate  under  scale 
mixtures  that  the  distribution  of  the  TSI.S  estimator  is  a  mixture  of  doubly 
noncentral  F  distributions  mixed  over  their  noneentra I i tv  parameters.  Fvi- 
tlenee  is  given  supporting  the  view  that  normal -theory  distributions  are 
robust  for  certain  mixtures.  In  Section  3  selected  percentiles  of  the  I'SI.S 
estimator  are  presented  for  contaminated  spherical  (Iaussian  and  spherical 
Cauchy  errors,  and  these  are  compared  with  the  standard  distributions. 

2.  NONSTANDARD  DISTRIBUTIONS 

bet  A(rx.s)  he  a  matrix  of  order  (r*s);  denote  by  L(Z)  =  N  (M,l  x£ ) 

r ,  s  -  r 

that  /.(r^s)  is  a  random  (Iaussian  matrix  with  expectation  F.(Z)  -  M,  the 
rows  of  7.  being  independent  with  the  same  dispersion  matrix  Z(s< s)  ;  and 

n 

let  (I  (M,l  xr)  he  a  scale  mixture  of  the  (Iaussian  laws  {N  (M,t""(I  xp)); 
r ,  s  -  r  -  r ,  s  -  -  r  - 

i  t.  (0,<«>)  1  whose  probability  density  function  (pdf)  has  the  representation 

g(Y)  =  (27i)‘rs/2^0°|T2rrr/2exp[-tr(Y-M)  '  (Y-M)  r_1 /2t  2  ]  d(l  ( t ) 

2 

with  d(r)  a  mixing  distribution  on  (0,°°) .  Further  let  \  (v,A)  he  the 


(2.  1) 


eh  i-squared  distribution  having  v  tie;’,  roes  of  I  ruudom  and  the  :  oticen  t  ra  I  i  ! 
parameter  A,  and  lot  I  ' (,  r ,  s  ,  A  ,  1 )  bo  tho  doubly  uonrenLral  i-  distribution 
with  numerator  parameters  (r,^)  and  denominator  parameters  (s,  ). 

Following  Anderson  and  Sawn  (1973),  consider  a  single  structural 
equation  of  the  tvpe 

V,  =  ry2  +  Vi  +  u  (;!- 

where  Vj  and  y.;  are  columns  of  the  observable  matrix  Y(T'2);  y  ^  (  |  •  K  )  is 
a  matrix  of  rank  K  T  cons  i  st  ing  of  known  exogenous  variable;;  ,  .tK.-ll 
is  .i  vector  of  parameters;  and  u(l'-l)  is  an  error  vector.  The  stru  turil 
equation  (2.2)  is  a  member  of  a  system  of  linear  equations  whose  reduced 
form  equations  are  represented  bv 

Y  =  Zli  +  V 

where  Z  =  (ZpZ,?)  is  a  matrix  of  order  (TxK)  and  rank  K  •  T  of  exi'twnoii,- 
variables;  1!  ( K  >;  2 )  is  a  matrix  of  reduced  form  parameters  partitioned  a 


conformably  with  Z;  and  V(Tx2)  is  a  matrix  of  random  d i st urbanees .  It  is 
.assumed  that  ’“22^  un^t  r;1nk  and  that  n0,;  /  0.  Anderson  and  Sawa 

(1973)  studied  the  k-class  estimators  p.  for  ;•  in  (2.2)  under  the  assumption 
that  L(V)  =  N  .,(0,  I,  x;>) ,  with  Si  =  a  (2-2)  matrix. 

1,2  1  -  i  I 

To  find  the  distribution  of  the  TSLS  estimator  for  B  under  scale  mix¬ 
tures  of  spherical  Gaussian  Jaws,  we  proceed  as  in  Anderson  and  Sawa  (1973, 
1979)  alter  cond i t ion  inn  on  the  mixing  parameter.  Their  reduction  to  a 
canonical  form  and  several  subsequent  transformations  are  valid,  condition¬ 
ally,  for  scale  mixtures.  Thus  starting  with  expression  (2.7)  of  Anderson 


and  Sawn  (1979),  wu  obtain  tin-  following  mixture  representation  for  t  he 
r  u"Hi  1 ,  1 1  i  vo  distribution'  f  unc  t  ion  (rdf)  ot  the  standard  i  red  I’SI.S  i-sl  iinator 


h 

conditioning  variable  ;  enter;;  tin.’  doubly  nmi'  cut  r;il  !•'  .1  i  ;  t  r  i  :  n  t  i  on  only 
through  its  noneent  ra  1  i  t  v  parameters.  I  he  uncond  i  t  i  ona  1  call  el  the 
standardized  ISIS  estimator  thus  is  a  mixture  ol  douhlv  noneent  ra.  1  !■'  dis¬ 
tributions  mixed  over  their  noneent ra 1 i t v  parameters,  i.e., 

2 ' *  ...  , 

l’(  -  — ~  (h-h)  ••  x)  =  /  l’(w.'w.  /w*w  •  it  (.  r )  !  t  i  dc  (  i  i  .  ( .' .  9 ) 

o  -  ()  *  I  ■  -  -  " 

It  is  clear  from  remarks  on  page  704  ol  Anderson  and  Sava  (197  1)  tint 
ptirtillel  results  ran  he  obtained  for  the  ordinary  least-squares  igii.S) 
estimator  b  .  Corresponding  to  (2.9)  wo  have 

!’(  “"(h  *  x)  =  /  1*  (w'w  /'w  'w  ■:  e  ( r )  |  t  )  dC  ( t  )  (2.10) 

1  -  o  ~  J~  *  •  -  -  - 

where 

l.( w'w  /w'w  )  =  I- (T-K,  ,T-K.  ,-..(  t  ),>..,( t  ))  .  (2.11) 

-  5-  j  •  4  4  1  1  1  2 

The  only  difference  between  (2.9)  and  (2.10)  is  the  number  of  decrees  of 
freedom  of  the  douhlv  noneent rn I  I’  distrihut  ion. 

We  conclude  tliis  section  with  a  brief  look  at  the  possible  robustness 
i > t  the  normal-theory  distribution  of  tile  TSI.S  est  imator  against  error  dis¬ 
tributions  arising  as  scale  mixtures  of  spherii-.il  Caussian  laws.  Anderson 

and  Sawa  (  1979)  provided  tables  for  the  edf  of  tho  standardized  TSI.S  eslinia- 

) 

tor  in  terms  of  the  degrees  of  freedom  (K.,)  and  the  parameters  t  and  .S“  de¬ 
fined  in  (2.7)  and  (2.4).  The  values  of  K.,  and  t  are  invariant  under  the 

>  °  2 

scale  change  u  •  ru  in  the  errors,  whereas  •  6  / t".  Thus  tiie  edf  of 
the  standardized  TSI.S  estimator  under  scale  mixtures  ot  Caussian  errors  is 
a  weighted  average  of  the  normal-theory  edf 's  corresponding  to  different 
values  of  5. 

A  perusal  of  the  tables  of  Anderson  and  Sawa  (1^79)  shows  that  tin- 


standardized  edf  often  remains  fairly  stable  over  a  wide  range  of  values  of 


/ 


Mixtures  over  these  rouges  accord  i  ng  1  y  .should  ■  j .  ■ !  i  i  .  •  1  . 

d  is  t  r  i  but  ions  for  which'  the  normal -t  heorv  d  i  st  rihut  ion  tin-  I'  !  ■  : 

is  robust  .  To  illustrate,  consider  the  two-point  mixing  bin:  i  sit  i  a 

■>  2 

i  (i“=l)  =  p  =  1  -  |'(i  “=5)  which  gives  a  contaminated  -.phii  i .  a  1  .an  ...  . 

error  distribution.  If  the  unconditional  1’Sl.S  estimator  ha.  it  -  .  . 

1 "  -  SOOO  with  a  fixed  parameter,  then  its  Pdf  is  a  mixture  . .  i  t;  n.  a  .  ! 

a  ■>  ;>  ’ 

theory  cdf's  cor  respond  in)’,  to  S"/\"  -  1000  when  ;  =  a  and  la  "  / .  '  iDtlo 

•> 

will'll  't  *"  =  1.  Selected  points  of  the  resultin'  edf  Were  evaluated  f  o  r 

t  •  {0,1, SI,  K.;  .  13,10],  and  p  ■' 0 .  a(»  ,0 . 7  5 , 0  . '»)  1  using  the  tallies  of 

Anderson  and  Sawa  (1970);  the  results  are  given  in  lahle  2.1  tor  K  t  =  3  and 

in  Table  2.2  for  K.,  -  it).  A  comparison  oi  the  norma  1 -theory  d  i  ,.l  r  i  but  i  on  t 

as  approximations  to  the  distributions  under  mixtures  confirms  the  robustness 
of  tin'  former  under  certain  types  of  mixtures. 

In  tile  following  section  we  undertake  a  broader  study  of  the  distribution 
of  tin'  standardized  TS1.S  estimator  under  various  error  distributions. 

3.  RKSUI.TS  AND  DISCISSION 

In  order  to  compute  the  edf  of  the  stand  i  rd i zed  ISI.S  estimator  under 
various  choices  for  the  mixing  distribution  (•(:),  we  follow  a  procedure  sim¬ 
ilar  to  the  one  sketched  al  the  end  of  the  preceding  section.  first,  lot- 
each  fixed  value  of  the  mixing  parameter  t ,  the  edf  of  tho  doubly  noucontral 
!■'  distribution  is  computed  using  the  algorithm  suggested  bv  Anderson  and 
Sawa  (1979).  The  resulting  values  of  these  rdf's  are  then  mixed  over  the 
relevant  values  of  r  according  to  ( f  (  t  )  .  Tin's  procedure  is  straightforward 
when  (1  ( t )  is  discrete.  for  continuous  mixtures  the  edf  C 1  ( t )  is  discretized 
at  the  second  step  on  partitioning  (0,«>)  into  intervals,  assigning  the  prob¬ 
abilities  to  mid-points  of  these  intervals,  and  proceeding  as  for  the  dis¬ 
crete  case.  This  procedure  is  expected  to  give  good  approximations  to  the 


ret  i  uni  su t  I  i  i  i  i  n t  I  v  . 


actual  distributions  when  tin1  discretization  is 

Our  numerical  studios  encompass  tliroo  di  nri-tf  mixtures  ol  Bo  i  non  i  i  i 
and  binomial  tvpo  giving  contaminat  od  spheric.il  danssi.m  error  ,  and  one 
continuous  mixture  using  an  inverse  o!ii  distribution  whiilt  -ives  ncri.nl 
Oaucliv  errors.  The  litter  is  an  example  ol  a  dist  r  i  but  ion  bavin-.',  .-xee  i  ve 

i 

tails.  In  particular,  our  Tv  no  1  eon  t  ami  na  L  i  on  uses  ‘i’ti'-li  0  .  v  > . 

■)  2  | 

!’(;"=  10)  =  0.05};  Type  2  eon  t  am  i  na  t  i  on  uses  |i'(i  -1)  “  0.5,  !'(  10) 

0.51;  and  Tvpo  3  coni  .iinin.it  inn  assigns  probabilities  to  i  .  1,2 . 16  i 

according  to  the  binonii.il  distribution 

0.05) '  ( 0 . a 5)  1  (  t.  !  ) 

with  r  =  i~  -  i.  These  disereto  mixtures  pose  no  difficulty  in  t  lie  rnmnu- 
ta  t  i  oils  . 

The  spherical  Cauchy  error  law  is  a  1-dimensional  spherical  Student'-, 
distribution  having  unit  degree  of  freedom;  multivariate'  Student's  errors 
were  considered  bv  Zell  nor  (1976)  in  connect imi  with  01.  S  estimation  in  a 

general  linear  model.  Spherical  Cauchy  errors  in  TS1.S  estimation  are 

-2  2  2 
obtained  from  L(t  ’)  =  y  (1,0)  on  evaluating  the  edf  of  \“(l,o)  over  226 

intervals,  assigning  the  probabilities  of  these  intervals  to  the  reciprocal 
of  their  midpoints,  and.  then  mixing  the  double  noncentr.il  K  distributions 
over  these  values.  As  doubling  the  number  of  intervals  bad  little  effect 
on  the  values  of  the  cdf  of  the  TS1.S  estimator,  the  choice  of  226  intervals 
tor  the  discrete  approximation  was  deemed  adequate.  The  tour  types  ol  mix¬ 
tures  studied  in  this  section  are  summarized  for  convenience  in  Table  3.1. 
Selected  percentiles  of  the  distributions  of  the  TSI.S  estimators  are 

given  in  Tables  3.2  -  3.7  along  with  the  interquartile  range  as  a  measure 

2 

of  concentration.  The  range  of  parameters  includes  6  t.  {100,300,1000, 
50001,  a  e  {0,1,5},  and  the  degrees  of  freedom  K  i 


{3,10}.  The  va  1  lies 


•rr-.irs  ar'c 


tabulated  by  Anderson  and  Sawa  (1979)  for  spin- ric.il  •'•ausaiaii  < 

given  along  with  our  own  eomput.ut  ions  for  t  hi-  l  hree  eontanin  at  ion  model;; 

and  for  spherical  Cauchy  errors.  The  medians  and  the  tail  probabilities 

reported  appear  to  he  accurate  within  +0.01,  while  the  interquartile 

ranges  mae  he  overstated  by  an  amount  in  the  range  ( 0 . 00 , 0 . 0 2 ) . 

Several  charac ter  is t i cs  of  the  1S1.S  estimator  were  noted  bv  Anderson 

and  Sawa  (1979)  in  the  ease  of  i  id  ''..mss  i. in  errors.  These  same  tendencies 

are  exhibited  under  the  additional  error  distributions  considered  here. 

hirst,  except  when  i  -  0,  the  TSI.S  estimator  has  negative  median  bias. 

This  hias  increases  with  i  and  with  the  degrees  of  iroedom,  K.)t  but  it 
•> 

decreases  with  "  .  Second,  tile  dispersion  of  the  TSI.S  estimator,  as 

measured  by  the  interquartile  range,  decreases  as  t  and  K.,  increase,  but 
a 

it  increases  witli  Together  these  two  observations  imply  that  as  t  or 

K  increase,  the  TSI.S  estimator  is  more  tightly  concentrated,  but  about 

values  further  removed  from  the  true  value  of  the  parameter.  The  same 

2 

holds  true  for  decreasing  values  of  ^  . 

in  the  preceding  section  we  noted  that  the  only  difference  between 
the  01, S  and  TSLS  estimators  lies  in  the  degrees  of  freedom  of  ihe  doubly 
noncentra]  F  distribution.  In  most  practical  circumstances  one  would 
expect  the  number  of  degrees  of  freedom  for  the  QI.S  estimator  to  be  greater 
than  that  for  the  TSI.S  estimator.  From  the.  foregoing  discussion  we  would 
then  expect  the  0LS  estimator  to  have  greater  bias  but  smaller  interquartile 
range  than  the  TSI.S  estimator.  Although  the  01. S  estimator  is  not  studied  hei 
numerically,  our  expectation  agrees  with  results  found  by  Raj  (1980)  in  a 
Monte  Carlo  study  using  independent  lognormal,  uniform,  and  double  expo¬ 
nential  errors. 

Comparing  properties  of  the  TSI.S  estimator  under  nonstandard  error 
distributions  to  those  for  Gaussian  errors  in  Tables  3.2  -  3.7,  we  note 


in 


that  the  median  bias  for  the  nonstandard  errors  is  never  less  than  that 
for  Caussian  errors.  Comparisons  of  the  interquartile  ranees  surest 
that  mixing  affects  location  of  the  derived  distribution  of  the  TSI.S 

estimator  to  a  greater  extent  than  scale. 

2 

For  large  values  of  the  parameter  it  is  seen  that  the  normal- 
theory  distribution  is  reasonably  robust  against  contaminated  Caussian 
errors  of  Types  1  and  3.  Type  2  contamination  and  Cauchy  errors  yield 
distributions  less  resembling  those  for  the  standard  case.  From  these 
studies  it  appears  that  rather  small  departures  from  the  standard  Caussian 
assumptions  are  not  crucial  over  a  modestly  wide  range  of  the  parameters 


of  the  distributions. 


of  coef f icients 


1  1 
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TABLK  2.1.  Distributions  of  the  standardized  TSLS  estimator  under  (a)  Gaussian  errors 
with  5‘/f  =  1000,  (b)  Gaussian  errors  with  6 '-  / 1 '  =  5000,  and  (<•;  ant uminatcd  Gaussian 
errors  obtained  as  a  mixture  of  (a)  and  (b)  with  proportions  ] —  p  and  p,  respectively, 
for  K.,  =  3  degrees  of  freedom. 


a  =  0 

a  =  1 

i  -  5 

X 

(a) 

(b) 

(c) 

p=0 . 50 

(a) 

(b) 

(e ) 

p=0 . 7  3 

(a) 

'b) 

(c  ) 

p=0 . 90 

-3.00 

0.0014 

0.0014 

0.0014 

0.0009 

0.0011 

0.0010 

0.0005 

0.0009 

0.0009 

-2.00 

0.0229 

0.0228 

0.0228 

0.0202 

0.0217 

0.0213 

0.0191 

0.0212 

0.0210 

-1  .00 

0.1587 

0.1587 

0.1587 

0.1664 

0.1612 

0.1620 

0.1667 

0.1623 

0. 1627 

-0.20 

0.4209 

0.4210 

0.4209 

0.4383 

0.4287 

0.4311 

0.4451 

0.4317 

0.4331 

0.00 

0.4997 

0.4997 

0.4997 

0.5175 

0.5077 

0.51 02 

0. 5245 

0.5109 

9.5123 

0.20 

0.5791 

0.5790 

0.5790 

0.5960 

0.5866 

0.5890 

0.6025 

0.5809 

0. 5910 

1.00 

0.8413 

0.8413 

0.8413 

0.8465 

0.8436 

0.8444 

0.8485 

0.8446 

0.6450 

2.00 

0.9771 

0.9772 

0.9772 

0.9750 

0.9762 

0.9759 

0.9742 

0.9758 

0.9757 

3.00 

0.9986 

0.9986 

0.9986 

0.9979 

0.9983 

0.9982 

0.9976 

0.9982 

0.9981 

TABLK  3.1.  Description  of  the  scale  mixtures  of  sphericaL  Gaussian  errors 
studied  numerically. 


Type 

Descr i p  t ion 

Type  1 

P(t2=I)  =  0.95 

;  I’(r2=10)  =  0.05 

Type  2 

P(t2=1)  =  0.5  ; 

l’(i2=10)  =  0.5 

Type  3 

P(x2=r+1)  =('r5) 

(0. 05) r (0. 95) 1 5-r 

r  e  {0,1,... 

.15} 

Hr  2)  =  x2O,0) 


Cauchy 


TABI.P.  3.2.  Selected  percer. t  i  1  es  and  range  of  the 
various  scale  mixtures  of  Gaussian  errors  for  K., 
various  A". 

1  TSI.S  estimator 
=  3  degrees  of 

under  Gaussian 
freedom,  a  =  0, 

and 

and 

I1RR0R 

DISTRIBUTIONS 

QUANTITY 

Gaussian 

Type  1 

Type  2 

Tv pe  3 

Cauchy 

6 2  =  100 

Med  ian 

0.00 

0.00 

0.00 

0.00 

0.00 

Interquart ile 

Range 

1.34 

1.34 

1.30 

1  .34 

1 .20 

2.5  Percentile 

-1.98 

CO 

CTv 

l 

-2.03 

-1.99 

-1.96 

97.5  Percentile 

1.98 

1.98 

2.03 

1  .99 

2.00 

o 

It 

Median 

0.00 

0.00 

0.00 

0.00 

0 . 00 

Interquart ile 

Range 

1.35 

1.35 

1  .  34 

1  .35 

1.26 

2.5  Percentile 

-1.97 

-1.97 

-1.98 

-1.97 

-1.96 

97.5  Percentile 

1.97 

1.97 

1.98 

1.97 

re 

O 

c 

5 2  =  1000 

Median 

0.00 

0.00 

0 . 00 

0.00 

0.00 

Intcrquarti le 

Range 

1 .  35 

1.35 

1  .  33 

I  .  33 

1  .  30 

2.5  Percentile 

-1.96 

-1.96 

-1.97 

-1.96 

-1.96 

97.5  Percentile 

1.96 

1  .  9t> 

1  .97 

1  .  96 

2 . 00 

=  5000 

Median 

0.00 

0.00 

0.00 

0.00 

0.00 

Interquart  ile 

Range 

1.35 

1  .  35 

1  .  35 

1  .  15 

1  .  34 

2.5  Percentile 

-1.96 

-1.96 

-1.96 

- 1  . 96 

-1.96 

97.5  Percentile 

1 . 96 

1.96 

1 .96 

1  .96 

1  .99 

TABLE  3.3.  Selected  percentiles  and  range  of  the  TSLS  estimator  under  Gaussian  and 
various  scale  mixtures  of  Gaussian  errors  for  =  3  degrees  of  freedom,  a  =  1,  and 
various  5^. 


ERROR 

DISTRIBUTIONS 

QUANTITY 

Gaussian 

Type  1 

Type  2 

Type  3 

Cauchy 

<52  « 

100 

Median 

-0.14 

-0.16 

-0.28 

-0.18 

-0.26 

Interquartile 

Range 

1.32 

1.34 

1.24 

1.33 

1.16 

2.5  Percentile 

-1 .84 

-1.83 

-1.80 

-1.81 

-1.82 

97.5  Percentile 

2.09 

2.08 

2.09 

2.11 

2.06 

62  = 

300 

Median 

-0.08 

-0.09 

-0.17 

-0.10 

-0.20 

Interquartile 

Range 

1.34 

1.35 

1.33 

1.34 

1.26 

2.5  Percentile 

-1.89 

-1.88 

-1.82 

-1.86 

-1.84 

97.5  Percentile 

2.04 

2.04 

2.08 

2.05 

2.06 

62  = 

1000 

Median 

-0.05 

-0.05 

-0.10 

-0.06 

-0.14 

Interquartile 

Range 

1.35 

1.36 

1.35 

1.35 

1.31 

2.5  Percentile 

-1.92 

-1.92 

-1.88 

-1.91 

-1.88 

97.5  Percentile 

2.00 

2.00 

2.04 

2.01 

2.04 

52  = 

5000 

Median 

-0.02 

-0.02 

-0.04 

-0.02 

-0.08 

Interquartile 

Range 

1.35 

1.36 

1.36 

1.35 

1.34 

2.5  Percentile 

-1.94 

-1.94 

-1.92 

-1.94 

-1.92 

97.5  Percentile 

1.98 

1.98 

2.00 

1.98 

2.02 

TABLE  3.4.  Selected  percentiles  and  range  of  the  TSLS  estimator  under  Gaussian  and 
various  scale  mixtures  of  Gaussian  errors  for  K.,  =  3  degrees  of  freedom,  a  =  5,  and 
various  6~ . 


ERROR 

DISTRIBUTIONS 

QUANTITY 

Gauss i an 

Type  1 

Type  2 

Type  3 

Cauchy 

6 2  =  100 

Median 

-0.19 

-0.22 

-0.40 

-0.2  5 

1 

o 

OJ 

CO 

Interquart i le 

Range 

1.30 

1.32 

1  .17 

1.31 

1 . 12 

2.5  Percentile 

-1.78 

-1.76 

-1.68 

-1.74 

-1 . 72 

97.5  Percentile 

2.13 

2.11 

2.07 

2.14 

2 . 04 

=  300 

Med  ian 

-0.  LI 

-0.12 

-0.2  3 

-0.  14 

-0.28 

Interquart i 1 e 

Range 

1.33 

1.34 

1  .  30 

1 .  34 

1.25 

2.5  Percentile 

-1.86 

-1.84 

-1.76 

-1.82 

-1  .  78 

97.5  Percentile 

2.06 

2.06 

2.10 

2.08 

2.06 

4  2  =  1000 

Med ian 

-0.06 

-0.07 

-0.13 

-0.08 

-0.20 

Interquart i le 

Range 

I  .33 

1.35 

1.35 

1  .  36 

1  .  32 

2 . 5  Percent i le 

-1.90 

-1.90 

-1  .83 

-1.88 

-1.84 

97.5  Percentile 

2.02 

2.02 

2.06 

2.03 

2.0b 

6?  =  5000 

Median 

-0.03 

-0.03 

-0.06 

-0.04 

-0.10 

Interquart ile 

Range 

1.35 

1  .35 

1  .  35 

1  .  36 

1.35 

2.5  Percentile 

-1.94 

-1.93 

-1.90 

-1.93 

-1.88 

97.5  Percentile 

1.99 

1.99 

2.01 

1  .99 

2.04 

TABLK  3.5.  Selected  percentiles  and  range  of  the  TSLS  estimator  under  Gaussian  and 

various  scale  mixtures  of  Gaussian  errors  for  K,  =  10  degrees  of  freedom,  a  -  0,  and 
r*>  J. 

various  o4- . 


QUANTITY 

Gaussian 

Type  1 

ERROR 

DISTRIBUTIONS 

Tv  pc  2 

1'ype  3 

Cauchy 

Med  ian 

0.00 

62  = 

0.00 

100 

0.00 

0.00 

0.00 

Interquartile 

Range 

L  .  29 

1.28 

1.13 

1  .  26 

1 .02 

2.5  Percentile 

-1.91 

-1.90 

-1  .75 

-1  .88 

-1.76 

97.5  Percent  i  Le 

1.91 

1.90 

1.75 

1.88 

1.78 

Median 

0.00 

6 2  = 

0.00 

300 

0 . 00 

0.00 

0.00 

Interquart il e 
Range 

1.33 

1.33 

1  .26 

1 .  32 

1.14 

2.5  Percentile 

-1.94 

-1.94 

-1  .88 

-1.93 

-1.84 

97.5  Percentile 

1.94 

1.94 

1 .88 

1  .93 

1  .87 

Median 

0.00 

32  = 

0.00 

1000 

0 . 00 

0 . 00 

0 . 00 

Interquart i le 
Range 

1.34 

1.34 

!  .32 

1.34 

1.24 

2.3  Percentile 

-1.96 

-1.95 

-1.93 

-1.95 

-1.90 

97 . 5  Percent ile 

1.96 

1.95 

1  .93 

1  .95 

1.92 

Med  ian 

0.00 

62  = 

0.00 

5000 

0.00 

0 . 00 

0.00 

Interquart i le 
Range 

1.35 

1.35 

1  .  35 

1.35 

1  .  30 

2.5  Percentile 

-1.96 

-1.96 

-1.96 

-1.96 

-1.93 

97.5  Percentile 

1.96 

1.96 

1.96 

1  .  96 

1.96 

TABLE  3.b.  Selected  percentiles  and  range  of  the  'l'Sl.S  estimator  under  Gaussian  and 
various  sjjale  mixtures  of  Gaussian  errors  for  K,  =  10  degrees  of  freedom,  u  =  1,  and 
various  <5". 


ERROR 

DISTRIBUTIONS 

QUANTITY 

Gauss  inn 

Type  ! 

Type  2 

Type  3 

Ca uchy 

S2  = 

100 

Median 

-0.59 

-0.62 

-0 . 89 

-0.72 

-0.71 

Interquart ile 

Range 

1.22 

1  .22 

1 .04 

1.17 

1.03 

2  .  5  i’ercent  i  le 

-2.19 

-2.19 

-2  .  20 

-2.21 

-2.14 

97.5  Percent  i  1 e 

1.44 

1.41 

1.10 

1.27 

1  .26 

■V  = 

300 

Med  ian 

1 

o 

LO 

O' 

-0.39 

-0  .  *)6 

-0.46 

1 

o 

00 

Interquart i Le 

Range 

1.30 

1.31 

1.23 

1  .29 

1.17 

2.5  Percentile 

-2.12 

-2.13 

-2.20 

-2.15 

-2.13 

97.5  Percentile 

1.69 

1.67 

1.42 

1  . 60 

1.52 

62  = 

1000 

Median 

-0.20 

-0.22 

-0.40 

-0.26 

-o .  42 

Interquart  i  le 

Range 

1  .  33 

1.34 

1  .32 

1  .  34 

1  .  30 

2.5  Percentile 

-2.06 

r-* 

O 

1 

-2.13 

-2.08 

-2.10 

9  7.5  Percent  ile 

1.83 

1.82 

1.67 

1 .  78 

1.71 

62  = 

5000 

Median 

-0.09 

-0. 10 

-0.  19 

-0.12 

-0.26 

Lnterquart  ile 

Range 

1.35 

1.35 

1  .  35 

1  .  35 

1.36 

2.5  Percentile 

-2.01 

-2.02 

-2.05 

-2 . 02 

-2.06 

97.5  Percentile 

1.91 

1.90 

1  .84 

1  .88 

1.85 

TABLE  3.7.  Selected  percentiles  and  range  of  the  TSLS  estimator  under  Gaussian  and 
various  scale  mixtures  of  Gaussian  errors  for  K.;  =  .10  degrees  of  freedom,  a  =  5,  and 
various  <5^. 


ERROR 

DISTRIBUTIONS 

QUANTITY 

Gauss i an 

Type  1 

Type  2 

Type  3 

Cauchy 

5 2  =  100 

Median 

-0.82 

-0.86 

-1 .28 

- 1  .  00 

-1.02 

Into i quart ile 

Range 

1.15 

1.16 

0.92 

1 . 08 

1.05 

2.5  Percentile 

-2.24 

-2.24 

-2.25 

-2.26 

-2 . 20 

97.5  Percentile 

1.19 

1.16 

0.82 

0.97 

1 . 04 

6-'  =  300 

Median 

-0.50 

-0.54 

-0.94 

-0.62 

-0.83 

Interquar t i le 

Range 

1.28 

l .  30 

1 .20 

1  .26 

1.22 

2.5  Percentile 

-2.17 

-2.18 

-2.24 

-2.20 

-2.18 

97.5  Percentile 

1.58 

1.55 

1.26 

1.44 

1.40 

5 2  =  1000 

Median 

-0.28 

-0.32 

-0.56 

-0.36 

-0.58 

Interquart ile 

Range 

1.33 

1.34 

1.32 

1  .32 

1.36 

2.5  Percentile 

-2.09 

-2.10 

-2.18 

-2.12 

-2.16 

97.5  Percentile 

1.77 

1.76 

1.56 

1.71 

1  .63 

6 2  =  5000 

Median 

-0.13 

-0.14 

-0.26 

-0 .  1  6 

-0.34 

Interquart i  le 

Range 

1.34 

1.35 

1  .  35 

1  .  34 

1.42 

2.5  Percentile 

-2.03 

-2.03 

-2.09 

-2.04 

-2.10 

97.5  Percentile 

1.88 

1.88 

1 .80 

1  .86 

1.81 
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distributions;  mixtures  of  distributions;  robustness. 

' 30  ABSTRACT  (Continue  on  revere*  elde  II  neceeaery  end  Identity  by  blink  number) 

(  in  estimating  the  coefficient  of  an  endogenous  variable  in  a  single 

equation  of  a  system  of  linear  equations,  Anderson  and  Sawa  (1971)  ex- 
■  pressed  the  distribution  of  the  two-stage  least-squares  (TSI.S)  estimator 

i 

j  as  a  doubly  noncentral  F  distribution.  We  relax  their  assumption  of  in- 
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20.  (Continued) 

dependent  Caussian  errors,  taking,  instead  a  scale  mixture  of  .spherical 
Caussian  laws  in  a  class  containing  the  spherical  stable  distributions. 
The  resulting  distribution  of  the  TS1.S  estimator  is  a  mixture  ol  doubly 
noncontral  F  distributions  mixed  over  the  noncentrality  parameters,  and 
lor  suitable  mixtures  the  normal-theory  distribution  is  robust.  Computa¬ 
tions  reported  tor  contaminated  spherical  Caussian  and  spherical  Cnuchv 


•rrors  are  compared  with  the  standard  case. 
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